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Background  Sky  Brightness  Measurements  for  Application  to  Space 
Surveillance  Systems 


Contract  FU9620-78-C-0013 

U.S.  Air  Force  Office  of  Scientific  Research 


Summary 


1.  Overall  objectives  and  strategy.  Studies  are  made  of  the  brightness,  color, 
and  angular  distribution  of  the  total  diffuse  astronomical  background  radiation 
and  of  its  components.  Ground-based  (lit.  Haleakala,  Hawaii),  Earth-orbital 
(Skylab),  and  deep-space  probe  (Pioneers  10  and  ll)  observations  of  zodiacal 
light  (multicolor,  brightness  and  direction  and  amount  of  polarization,  different 
times  of  year  and  different  locations  in  the  solar  system,  all-sky  coverage  to 
within  30®  of  the  sun)  are  used  to  derive  information  on  the  physical  and 
optical  properties  of  interplanetary  dust,  its  origin  and  dynamics,  its  location 
in  space,  and  its  contribution  to  the  infrared  astronomical  background.  In  a 
search  for  interrelationships,  results  are  compared  with  meteoritic  data,  with 
photometric  data  on  comets,  and  with  properties  of  interstellar  grains. 

2.  Accomplishments  to  date.  Background  starlight  has  been  mapped  in  blue  and  red 
wavelengths  over  the  entire  sky,  using  Pioneer  10  observations  from  beyond  the 
asteroid  belt  where  the  zodiacal  light  is  negligible.  Zodiacal  light  has  been 
isolated  in  observations  between  1  and  3  AU  to  derive  heliocentric  dependence. 

The  general  solution  for  inversion  of  the  zodiacal  light  brightness  integral 
has  been  derived  and  is  used  to  yield  the  scattering  function  of  interplanetary 
dust  for  each  position  traversed  by  the  probe  moving  through  the  solar  system. 
Properties  of  the  dust  are  found  to  change  with  heliocentric  distance.  Polari¬ 
zation  of  zodiacal  light  was  determined  at  the  north  celestial,  south  ecliptic, 
and  north  galactic  poles.  For  inner  zodiacal  light,  symmetry  plane  of  inter¬ 
planetary  dust  is  close  to  orbital  plane  of  Venus,  At  110°  from  the  sun,  it 

is  close  to  orbital  plane  of  Mars.  There  appears  to  be  more  than  one  symmetary 
"plane"  for  the  dust;  position  of  maximum  dust  density  relates  closely  to  orbital 
planes  of  the  planets.  Mini-workshop  held  in  Albany  to  initiate  Joint  studies 
of  zodiacal  light  and  comet  data  and  interstellar  grains.  Dust  rotation/bursting 
studied  as  possible  important  mechanism  in  comets,  interplanetary  space,  reflection 
nebulae,  stellar  envelopes,  interstellar  medium.  Code  completed  for  complete 
solution  to  electromagnetic  scattering  by  spheroids.  Started  comparison  of  ob¬ 
servations  of  Comet  1965  VIII  with  scattering  measurements  of  irregular  particles 
and  with  model  calculations  for  spheres,  spheroids,  and  cylinders  to  derive 
spatial  distribution  and  physical  properties  of  dust  in  tail.  Zodiacal  light 
polarization  observations  on  antisolar  hemisphere  contest  claim  by  others  that 
interplanetary  dust  is  aligned  and  that  alignment  changes  over  short  times.  In- 
house  facility  completed  for  digitizing  entire  library  of  ground-based  observations 
from  Mt.  Haleakala,  I96U-I968.  Night  sky  observing  facility  at  Mt.  Haleakala 


1 


v> 


m 


Background  Information 

In  October  1976  we  began  an  AFOSR- funded  program  to  complete  the  re¬ 
duction  and  analysis  of  photometric  data  from  the  Pioneer  10/11  deep-space 
probes  and  from  Skylab  and  to  combine  these  data  with  ground-based  results 
and  theoretical  studies  to  map  the  diffuse  astronomical  background  radiation  • 
background  starlight  and  zodiacal  light  -  over  the  sky  and  to  determine  the 
nature  of  the  interplanetary  dust.  This  report  summarizes  results  obtained 
through  the  third  year  of  this  program  and  outlines  ongoing  tasks  in  the 
fourth  year,  which  is  supported  by  Grant  AFOSR- 8O-OOU3.  Additional  information 
is  given  in  the  March  1978  Final  Report  on  Grant  AFOSR-76-3091  and  in  the 
March  1979  Interim  Scientific  Report  on  the  subject  contract. 

Summary  of  Program  Objectives,  FY  1978-1979 


1.  To  complete  the  all-sky  mapping  of  Pioneer  10/11  data  on  background  star¬ 
light  in  blue  and  red  wavelengths,  using  observations  from  beyond  the 
asteroid  belt  where  the  zodiacal  light  is  negligible. 

2.  To  define  the  brightness  of  zodiacal  light  as  a  function  of  heliocentric 
distance,  elongation  (angular  distance  from  the  sun),  and  ecliptic  lati¬ 
tude. 

3.  To  determine  the  polarization  of  zodiacal  light  in  the  ecliptic  and  at 
the  north  and  south  celestial,  ecliptic,  and  galactic  poles  -  as  a 
function  of  heliocentric  distance. 

4.  To  establish  as  far  as  possible  the  nature  of  the  particles  giving  rise 
to  the  zodiacal  light. 

5.  To  continue  the  study  of  electromagnetic  scattering  by  spheroids  and 
cylinders  and  to  compare  these  results  to  Mie  theory. 


In  subsequent  sections  we  outline  the  status  of  these  objectives  and  of  related 
topics  as  of  the  end  of  FY  1979-  We  also  list  papers  which  were  published  or 
are  in  preparation  based  on  AFOSR  support. 


Pioneer  10/11  Results 

1.  Background  Starlight.  Between  March  1972  and  February  1976,  76  full  or 
partial  sky  maps  were  made  by  Pioneer  10.  Pioneer  11,  starting  in  April  1973, 
made  67  sky  maps.  Pioneer  11  was  retargeted  to  encounter  Saturn  in  1979; 
Pioneer  10  is  continuing  on  a  path  that  will  make  it  the  first  probe  to  leave 
the  solar  system.  Sky  mapping  observations  from  both  spacecraft  ceased  in 
May  1976.  For  measurements  made  between  1  and  2.8  AU,  sky  maps  contain  data 
on  background  starlight  and  zodiacal  light.  Since  the  zodiacal  light  was 
found  to  be  negligible  beyond  2.8  AU,  these  sky  maps  contain  data  on  back¬ 
ground  starlight  alone.  Our  procedure  was  to  use  merged  Pioneer  10/11  back¬ 
ground  starlight  data  from  a  number  of  maps  beyond  2.8  AU  to  isolate  the 
zodiacal  light  in  maps  at  smaller  heliocentric  distances.  98$  of  the  sky  is 
now  "mapped"  in  blue  and  red  wavelengths  in  a  brightness  contour  format 
similar  to  that  shown  in  the  March  1978  Final  Report  on  FY  1978  activities. 
The  remaining  sky  regions  have  been  processed  through  all  but  the  last  step 
prior  to  creating  final  contour  plots. 

As  noted  in  earlier  reports,  high  resolution  maps  of  the  background  star¬ 
light  in  several  colors  are  needed  to  separate  the  major  components  of  the 
light  of  the  night  sky  in  observations  made  from  within  the  zodiacal  cloud 
and  are  also  important  from  the  point  of  view  of  knowing  the  overall  topology 
of  the  diffuse  astronomical  background  radiation.  These  maps  are  also  the 
starting  point  in  the  separation  and  analysis  of  the  background  starlight 


components  themselves  (integrated  starlight,  diffuse  galactic  light,  cosmic 
light) . 


In  response  to  a  number  of  requests  for  these  data,  we  have  assembled 
a  background  starlight  data  file.  This  file  or  data  set  is  composed  of  3 
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parts: 

(1)  A  tabulated  matrix  in  galactic  coordinates  (b,l)  of  background 
starlight  in  two  colors  over  the  sky,  in  2-deg  increments  of  b,l. 

(2)  Contour  plots  of  blue  and  red  brightnesses  to  the  scale  of  Becvar's 
Atlas  Coeli.  A  sample  of  these  data  is  scheduled  for  publication 
in  the  May  1980  issue  of  Sky  and  Telescope . 

(3)  "Posted"  numbers  for  each  field  of  view  data  number  to  construct 
the  contours  in  (2)  above.  This  allows  one  to  determine  the 
orientation  of  the  Pioneer  field  of  view  for  each  data  number 
over  the  sky. 

With  these  parts,  a  user  could  determine  what  his  instrument  would  see  in 
absolute  units  for  any  size,  shape,  and  field  of  view  orientation  -  to  the 
Pioneer  10/11  instrument  limit  of  V  magnitude“6,  which  is  typical  of  most 
large  field  photometric  studies  of  the  diffuse  astronomical  background 
radiation.  These  data  are  already  being  provided  on  a  limited  basis  to  the 
Max-Planck-Institut  flir  Astronomie,  Heidelberg,  and  to  the  Tokyo  Astronomical 
Observatory.  Sample  data  files  were  prepared  for  display  and  explanation  at 
the  Montreal  General  Assembly  of  the  IAU  -  for  the  report  of  Commission  21' s 
Background  Starlight  Working  Group  and  for  IAU  Commissions  28  (Galaxies),  33 
(Structure  and  Dynamics  of  the  Galactic  System),  3^  (interstellar  Matter  and 
Planetary  Nebulae) .  Based  on  the  large  number  of  requests  for  these  data, 
we  are  investigating  costs  associated  with  having  these  data  commercially 
reproduced. 

2.  Zodiacal  Light.  All  pre-Pioneer  10/11  observations  of  the  zodiacal 
light  yield  only  the  brightness  due  to  the  scattering  of  sunlight  by  a 
column  density  of  interplanetary  dust.  Observations  from  the  Pioneer  10/11 
deep-space  probes  were  unique  in  that  the  spacecraft  could  be  used  as  optical 
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probes.  The  resulting  zodiacal  light  observations  made  in  the  plane  of  the 
spacecraft  trajectory  (approximately,  but  not  exactly,  the  ecliptic)  could 
be  inverted  to  yield  the  brightness  per  unit  line  of  sight  as  a  function  of 
scattering  angle.  This  exact  method  of  inversion  (Dumont,  Planet .  Space  Sci., 
21,  2149,  1973)  was  first  applied  to  Pioneer  10  observations  as  described  by 
Schuerman  ( Space  Research  XIX,  447-450,  1979)*  However,  the  trajectories  of 
Pioneer  10/11  were  not  exactly  planar  as  required  by  Dumont's  2-dimensional 
inversion  formula;  the  inclinations  of  the  spacecraft  orbits  varied  up  to 
15°  with  respect  to  the  ecliptic.  We  have  thus  generalized  the  inversion 
formula  to  three  dimensions  (Schuerman,  Planet .  Space  Sci. ,  27,  551-556,  1979, 
attached),  and  this  3-dimensional  inversion  methodology  has  been  applied  to 
Pioneer  10  data  at  19  heliocentric  distances  ranging  from  1  to  3  AU.  Final 
results  were  presented  at  IAU  Symposium  90,  Solid  Particles  in  the  Solar 
System,  Ottawa,  August  1979  (abstract  attached).  Of  particular  interest  is 
the  result  that  observations  from  Pioneer  10  do  not  verify  the  two  major  long¬ 
standing  assumptions:  (l)  the  dust  density  goes  as  some  power  of  the  helio¬ 
centric  distance,  and  (2)  the  dust  properties  are  independent  of  location. 

An  interesting  by-product  has  emerged  from  our  consideration  of  the 
symmetry  of  the  interplanetary  dust  complex.  The  inversion  method  mentioned 
above  relies  on  the  assumption  of  cylindrical  symmetary  of  the  dust  cloud. 

In  looking  for  mechanisms  which  might  violate  this  symmetry,  we  have  derived 
the  general  solution  for  the  restricted  3-body  problem  which  includes  the 
effect  of  radiation  pressure  on  the  third  (negligible  mass)  body.  The  positions 
of  the  equilibrium  ( Lagrangian)  points  are  dramatically  altered  by  radiation 
pressure;  no  overabundance  of  Bmall  particles  is  now  expected  to  be  trapped 
at  the  L4  and  L5  points  of  sun-planet  systems.  Details  of  the  3-tody  problem 


with  radiation  pressure  were  presented  at  Ottawa  and  are  in  press  (abstract 
attached);  certain  solutions  may  be  of  importance  for  both  solar  system  dust 
dynamics  and  interstellar  star  formation. 

Related  Results 
-  Skylab. 

The  availability  of  Pioneer  10/11  all- sky  data  on  background  starlight 
makes  it  possible  to  accurately  isolate  zodiacal  light  in  our  Skylab  and 
ground-based  observations.  The  use  of  this  method  and  these  data  has  re¬ 
sulted  in  consistently  good  intercomparisons  among  our  own  results  and  with 
the  results  of  Dumont  and  Sanchez.  This  is  especially  noteworthy  in  that 
the  Skylab,  Pioneer  10/11,  ground,  and  Dumont  and  Sanchez  measurements  were 
made  with  completely  different  instruments,  using  different  techniques  of 
measurement  and  data  reduction  and  completely  independent  absolute  cali¬ 
brations  . 

In  this  manner,  the  total  and  polarized  brightnesses  and  the  degree  of 
polarization  of  zodiacal  light  were  determined  for  all  fixed  position  Skylab 
observations:  at  the  north  celestial,  south  ecliptic,  and  north  galactic 
poles  and  at  90°  from  the  sun  in  the  ecliptic  (vernal  equinox);  the  results 
were  presented  at  Ottawa  (see  attached  abstract).  The  same  methodology  is 
now  being  applied  to  all  of  the  sky-scanning  observations  obtained  from 
Skylab.  The  polarization  data  are  reduced  in  a  different  manner;  they  are 
not  affected  by  background  starlight  except  at  low  galactic  latitudes. 
Polarization  results  to  date  confirm  our  earlier  finding  that  the  directions 
of  polarization  of  zodiacal  light  are  symmetric  with  respect  to  the  sun/ 
anti  sun,  with  the  result  that  the  particles  do  not  have  to  be  aligned. 
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-  Ground-based  data. 

During  the  report  period  existing  ground  observations  continued  to  be 
used  to  provide  information  on  the  symmetry  plane  of  the  zodiacal  cloud. 

Our  results  have  shown  that  the  symmetry  plane  is  near  the  orbital  plane 
of  Venus  at  elongations  from  32  to  50°  from  the  sun  (Misconi  and  Weinberg, 

Science,  200,  1484-1485,  1978).  Sky lab  data  at  110  to  l40°  from  the  sun 
showed  that  the  symmetry  plane  is  between  the  orbital  plane  of  Mars  and  the 
ecliptic  plane  at  those  elongations  (Misconi,  et  al..  Bull.  MS,  9,  620,  1977)- 
These  results  lead  us  to  believe  that  there  is  no  unique  plane  of  symmetry 
for  the  zodiacal  cloud,  but  rather  a  "multiplicity  of  planes"  associated  with 
the  orbital  planes  of  the  planets.  Additional  observational  data  are  needed, 
especially  at  elongations  beyond  50°.  For  this  purpose  we  are  reactivating 
our  night-sky  facility  atop  10,000  ft  Mt.  Haleakala,  Maui,  Hawaii.  The  instrument 
that  will  be  used  is  a  l4-color  (near  UV  to  near  IR)  photoelectric  polarimeter 
with  its  own  alt-azimuth  mounting  and  control  and  data  recording  system. 

-  Digitization  of  existing  ground-based  observations. 

In  the  period  1965  to  1968  we  made  photopolarimetric  observations  on 
more  than  450  nights.  The  major  portion  of  these  data,  and  all  single-color 
(5300A)  observations  dating  back  to  1961,  were  recorded  in  analog  (strip-chart) 
form;  funds  were  not  available  at  that  time  to  acquire  and  perfect  a  digital 
recording  system.  Approximately  20$  of  these  strip-chart  recordings  have 
been  digitized  using  manual  techniques,  a  semi-automatic  analog- to-digital 
chart  reader,  and  a  commercially- available  digitizing  service  which  utilizes 
a  high  speed,  laser  scanning  digitizer  (i/o  Metrics,  Palo  Alto,  Calif.). 

These  methods  are  now  prohibitive  in  cost,  and  we  have  developed  an  in-house 
capability  for  accurately  digitizing  these  strip-chart  recordings.  Several 
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years  ago  the  University  acquired  a  computer- controlled  TV  camera  (EMR  Optical 
Data  Digitizer)  as  the  principal  ingredient  of  a  digitizing  service  for  Uni¬ 
versity  and  local  research  groups.  We  have  refurbished  and  adapted  this  digiti¬ 
zer  for  use  with  photographic  negatives  of  our  strip-chart  recordings;  the 
system  is  now  operational,  reproducible,  and  accurate.  We  are  in  the  process 
of  photographing  those  recordings  bearing  directly  on  the  photometric  axis 
of  zodiacal  light  and  on  other  aspects  of  our  current  studies. 

-  Electromagnetic  scattering  by  spheroids. 

The  complete  solution  for  the  scattering  of  plane  electromagnetic  waves 
incident  upon  homogeneous  spheroids  was  given  by  Asano  and  Yamamoto  (Applied 
Optics ,  14,  29,  1975)-  These  first  solutions  contained  lengthy  and  involved 
mathematics  and  calculations  and  were  illustrated  only  for  the  case  of  a 
purely  real  refractive  index.  R.  Schaefer  has  worked  through  the  entire 
mathematics  associated  with  the  spheroid,  as  part  of  his  Ph.D.  thesis  re¬ 
search  at  SUNY  Albany,  and  he  has  extended  their  original  work  to  include  both 
dielectric  and  absorbing  particles.  The  computer  code  was  completed  and  docu¬ 
mented  during  the  report  period,  and  a  number  of  particle  and  ensemble  calcu¬ 
lations  have  been  made.  In  studies  just  completed,  excellent  agreement  was 
obtained  in  the  first  comparison  of  spheroid  calculations  with  measurements 
in  our  microwave  analog  scattering  laboratory. 

Additional  topics  in  process  include  dust  rotation/bursting,  development 
of  joint  (Alb any- Leiden- Seoul)  studies  of  zodiacal  light  and  comet  data  and 
interstellar  grains,  the  search  for  color  effects  in  the  brightness  and 
polarization  of  zodiacal  light,  an  observational  model  of  the  inner  zodiacal 
light,  the  reduction  of  Palomar  Sky  Survey  star  count  data  and  their  com¬ 
parison  with  Pioneer  10/11  background  starlight,  and  models  of  the  diffuse 
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Optical  detection  of  local 
,  •  j  interstellar  particles 

.i  during  an  out-of-ecliptic  mission 

The  contribution  to  the  zodiacal  light  (ZL)  by  interstellar  grains 
has  been  estimated1- *,  but,  as  suggested  here,  new  observations 
and  analyses  of  the  ZL  and  of  the  local  interstellar  medium 
provide  a  firmer  basis  for  placing  limits  on  the  brightness  of  the 
interstellar  (IS)  relative  to  the  interplanetary  (IP)  component. 
The  proposed  NASA/ESA  ‘out-of-ecliptic’  (or  ‘solar  polar') 
mission4  provides  the  best  and  perhaps  the  only  observational 
method  to  distinguish  between  these  two  components. 

To  estimate  the  interstellar  dust  density,  we  make  the  usual 
assumption  that  the  number  density  of  the  interstellar  particles 
(/in)  is  related  to  the  local  hydrogen  number  density  in  the  ratio 
of  10~n  (see  ref.  5).  The  local  hydrogen  density  within  a  few 
parsecs  of  the  Sun  has  been  observed  to  be  ~0.1  cm-3  (ref.  6) 
which  suggests  that  we  are  currently  passing  through  an  inter¬ 
cloud  region  of  the  galaxy.  We  infer  from  this  gas  density  that  the 
IS  dust  content  is  about  %» I0“13  cm"1. 

The  effect  of  the  Solar  System  on  the  ambient  IS  complex  of 
gas  and  dust  is  unknown.  Although  several  mechanisms  have 
been  proposed  either  for  capture7,4  or  the  sweeping  out  of 
interstellar  grains,  none  of  the  results  are  conclusive  because  of 
uncertainties  in  the  basic  assumptions.  The  most  efficient 
mechanism  for  deflecting  interstellar  particles  (if  charged)  is 
probably  the  Lorentz  force  exerted  by  the  magnetic  field  asso¬ 
ciated  with  the  solar  wind9.  Levy  and  Jokipii10  estimate  the  ratio 
of  the  Lorentz  to  the  gravitational  force  ( FJF0 )  at  1  AU  to  be 
10,  but  their  use  of  an  IS  particle  size  of  0.02  |un  is  too  small  by 
at  least  a  factor  of  five  for  the  classical  size  particles".  If  a 
realistic  particle  size  of  characteristic  radius  0.12  pm  is  used,  a 
value  of  F,/Fa  of  the  order  of  0.4  (rather  than  >1)  is  derived, 
and  sweeping  out  of  the  IS  particles  should  not  occur.  Given  the 
absence  of  detailed  dust  orbit  calculations  as  well  as  the 
uncertainties  in  both  the  initial  grain  charge  properties  and  in 
the  possibility  of  the  grains  discharging  in  the  IP  medium,  we 
base  our  calculations  on  the  assumption  that  the  grains  are 
neither  accreted  nor  swept  out,  but  that  the  IS  grains  in  the  Solar 
System  are  uniformly  distributed  with  a  density  characteristic  of 
the  local  IS  medium. 

We  estimate  the  monochromatic  IS  contribution  to  the  ZLby . 
scattering  of  solar  radiation  from  an  arbitrary  distribution  of 
interstellar  particles  from  the  brightness  integral11 


which  becomes,  for  a  uniform  IS  density, 
t  _  _  /Dif  . 


Zu^  ttafaFo 
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where  LOS  is  line-of-sight,  /o( R&/R)*  is  the  solar  flux  on  a 
typical  volume  element  at  a  heliocentric  distance  R,  0  is  the 
scattering  angle,  ara(0)  is  the  differential  scattering  cross  section 
of  the  IS  dust,  and  df  is  measured  along  the  LOS.  The  integral  in 
equation  (lb)  is  calculated  from 


where  b  is  the  closest  solar  approach  of  the  LOS,  c  is  the 
elongation  angle  (see  Fig.  1),  and  k  is  the  wave  number  (2tt/A) 
of  the  incident  radiation.  For  the  scattering  function  Fffl),  we 
use  a  dielectric  dust  model”.  At  a  wavelength  of  4,400  A,  the 
conversion  of  cgs  radiance  units  to  the  unit  Sl0(V)  (using  the 
observational  parameters  recommended  by  Sparrow  and 
Weinberg14)  is 


s,o(vy 


•  1.2x  1 0*9  erg  cm-*  s" 


where  S|0(V)  is  defined  as  the  radiance  produced  by  a  10th 
visual  magnitude  solar  (G2V)  star  per  square  degree. 

According  to  the  model  shown  in  Fig.  1,  at  the  Earth  (b  — 

1  AU  and  fore  «■  90“X  the  IS  contribution  to  the  ZL  is  only  about 
1  S.o(VO-  This  value  is  not  only  below  the  limit  of  detectability  of 
~2-3S»(V)  (ref.  IS);  it  is  also  totally  masked  by  the 
200  Sto(V)  contributed  by  the  IP  particles.  On  the  other  hand, 
viewing  sunwards  from  a  space  craft  (located  at  R  **  1  AU)  along 
an  LOS  parallel  to  the  ecliptic  plane  and  1/2  AU  above  it 
(out-of-the-edipticX  the  IS  brightness  estimated  from  equation 
(2)  for  the  then  appropriate  value  of  e  *»  27*  is  —28  Sio(  10  which 
is  well  above  the  limit  of  detectability.  Furthermore,  this  IS 
contribution  may  be  distinguishable  from  the  IP  component. 

An  observational  upper  limit  to  the  IS  brightness  component 
is  established  by  noting  that  the  ZL  is  close  to  solar  colour  from 
about  2,200  A  to  2.4  pm  (ref.  16  and  refs  thereinX  Since  the  IS 
particle  scattering  is  colour  (wavelength)  dependent  in  this 
range,  this  may  be  used  to  establish  an  upper  limit  to  the  IS 
component.  Our  calculations  are  confined  to  the  wavelength 
spread  of  the  UBV  system  where  more  data  are  available  for 
separating  the  ZL  and  the  background  starlight.  The  ZL  in  this ' 
wavelength  region  is  known  to  be  solar  colour  to  within  about 
5%. 

Assume  that  the  ZL  consists  of  two  components,  the  IP  part 
which  is  solar  colour  and  the  IS  component  whose  wavelength 
dependence  follows  the  IS  dust  model  previously  cited.  Thus  if 
we  consider  the  two  wavelengths  5,400  A  and  3,600  A,  the  IP 
component  has  a  spectral  ratio  /ip(5,400)//ip(3,600) = 
1(/]P^/ip(A);  I  is  normalised  to  solar  colour).  According  to 
Greenberg  and  Hanner31,  the  IS  component  has  a  spectral  ratio 
7B(5,400)//,s(3,600)= 0.77  at  e  =  90*. 

The  5%  colour  variation  constraint  then  implies 

/iP+/,s(3.600)  1  +/is(3.600)//|p  .. 

/,p+  /is(5,400)  "  1  +  0.77  x  /ls(3,600)//,p  1  J 


so  that 


75(3,600)«0.25  or  -^«0.20 

/IP  izL 


This  means  that  no  more  than  about  20%  of  the  total  zodiacal 
light  as  viewed  in  the  ecliptic  and  at  90°  elongation  can  be 
produced  by  scattering  off  of  interstellar  type  particles. 

Fig.  1  Resulting  ZL  brightness  in  St0  (V)  units  per  unit  (au) 
distance  of  closest  solar  approach  along  the  line  of  sight  from  )-< 
to  $  =  it.  Inset  shows  appropriate  geometry. 
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As  knowledge  of  the  size  and  spatial  distribution  of  the  IP 
particles  out  of  the  ecliptic  is  lacking,  we  derive  the  maximum 
amount  of  IP  dust  J  AU  off  the  ecliptic  that  would  still  allow  the 
IS  component  to  be  distinguished.  Let  f-nlr  be  this  amount, 
where  n,r  is  the  dust  density  in  the  neighbourhood  of  the  Earth 
and /  is  the  fraction  to  be  determined.  As  a  reference  point,  the 
neighbourhood  of  the  Earth  has  a  volumetric  scattering  cross 
section  of  <rtr( 90*)  =  7.8  X 10'”  cm"*  sr-1  (ref.  17).  For  0  - 
45°,  the  scattering  cross  section  o(r(4S0)  is  about  1.77  o-ir(90a) 
(ref.  18). 

Our  minimum  distinguishability  criterion  is  that  the  ratio  of 
the  interstellar  to  interplanetary  contributions  be  greater  than 
25% .  For  the  case  treated  above  for  the  interstellar  contribution 
at }  AU  above  the  ecliptic,  then 

Zn  ms<r»(4S*)  7.5xlQ-'lXlQ-” 
Z0.’"ntra,W)mfX  1.77x7.8x10'” 

from  which  we  obtain /<  0.20. 

This  value  of  /  is  close  to  but  smaller  than  those  derived  from  a 
series  of  models  by  Giese**  which  produced  results  consistent 
with  those  inferred  from  observations  from  the  Earth.  Most  of 
these  models  predict  the  dust  density  (assumed  to  be  purely 
interplanetary)  to  decrease  to  20%  between  heights  of  0.5  and 
1  AU  above  (or  below)  the  ecliptic. 

These  rough  calculations  suggest  there  is  a  real  possibility  of 
optically  detecting  a  local  IS  dust  component  in  the  Solar 
System.  Although  our  criterion  of  distinguishability  is  based  on 
colour  arguments  alone,  both  polarisation  and  scattering  angle 
dependence  provide  additional  tools  to  distinguish  between  the 
larger  (IP)  and  smaller  (IS)  particles.  The  chance  of  isolating  the 
IS  component  increases  with  decreasing  elongation  angle. 
Observations  shoud  be  made  along  a  line  of  sight  which  is 
parallel  to  the  ecliptic  arid  is  as  close  to  the  Sun  as  possible. 
Fortunately,  such  observations  can  also  be  inverted;  not  only 
can  the  total  brightness  ( Z )  be  obtained,  but  when  the  spacecraft 
is  near  maximum  height  the  value  no-  can  be  directly  determined 
from  the  measurements10. 

Assuming  that  a  positive  identification  of  interstellr  dust  will 
be  made,  we  are  on  the  threshold  of  observing  physical 
-  parameters  of  these  particles  which  could  previously  only  be 
inferred  from  observations^!  stellar  extinction  and  scattering  in 
reflection  nebula,  which  are  normally  hundreds  of  parsecs  from 
our  local  neighbourhood.  For  example,  given  a  slightly  more 
optimistic  estimate  (by  perhaps  a  factor  of  2)  of  the  local 
interstellar  dust  density,  we  might  be  able  to  observe  how  they 
are  affected  by  interaction  within  the  Solar  System. 
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Extent  of  hot  intergalactic 
gas  In  the  duster  Abell  2218 

A  SMau.  cooling  of  the  microwave  background  towards  the 
centre  of  the  rich,  distant  cluster  of  galaxies  Abell  22 18  has  been 
detected1*1.  This  decrement  has  been  attributed  to  the  inverse 
Compton  scattering  of  photons  of  the  background  radiation  by 
electrons  in  a  hot  intracluster  plasma,  as  predicted  by  Sunyaev 
and  Zel’dovich4.  We  report  here  new  observations  in  several 
different  directions  towards  the  duster  which  confirm  the  earlier 
results  and  provide  an  estimate  of  the  extent  of  the  gas. 

The  new  measurements  were  made  during  the  period  Febru 
ary  1977  to  May  1978  at  the  SRC  Chilbolton  Observatory,  and 
consist  of  199  h  of  observation  at  points  north  and  south  of  the 
Abell  cluster  centre1  and  54  h  further  at  that  position.  As  in 
previous  measurements*'1,  these  observations  were  made  at 
prime  focus  of  the  25-m  telescope  operated  at  10.6  GHz,  using  a 
twin-beam  system  with  a  beam-throw  of  15  arc  min  In  azimuth 
and  HPBW  of  4.5  arc  min.  The  receiver  incorporated  a 
parametric  amplifier  with  noise  temperature  about  180K  and 
bandwidth  250  MHz.  The  experimental  method  and  data 
reduction  have  been  described  in  detail  elsewhere*. 

Table  .1  shows  the  coordinates  of  the  positions  observed  and 
the  temperature  decrements  measured  there.  No  corrections  for 
the  presence  of  nearby  radio  sources  are  needed0.  The  present, 
results  are  consistent  with  earlier  measurements  near  dec  66* 
20' — as  discussed  by  Birkinshaw.  Gull  and  Northover*.  The 
measurement  at  dec  66*  18’  was  made  to  check  on  the  strong 
detections  at  66*  16*  and  66°  20'  and  to  furnish  another  point 
close  to  the  cluster  centre.  Its  presence  leads  to  an  appreciable 


Fig.  1  The  galaxy  strip  count.  N,  in  arc  min  bins  and  the  observed 
temperature  decrement  at  each  point,  A T,  plotted  as  functions  of 
declination .  S  (1 950.0 coordinates).  The  solid  line  is  the  AT  profile 
of  the  best-fitting  y  •»  5/3.  ( =»  -0.15  atmosphere  convolved  with 
the  beamshape  of  the  Chilbolton  telescope.  The  linear  scale,  D,  at 
the  distance  of  the  cluster  is  also  marked. 
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INVERTING  THE  ZODIACAL  LIGHT  BRIGHTNESS  INTEGRAL 

(Received  11  September  1978) 

Ah>nl — Considerations  of  the  geometry  appropriate  to  obaervatiom  of  the  zodiacal  light  made  from 
out  of  the  ediptic  plane  yield  the  general  inversion  of  the  brightness  integral.  The  brightness  per  unit 
volume  of  interplanetary  space  can  thus  be  determined  in  the  immediate  neighborhood  of  the 
spacecraft  in  directions  confined  to  a  unique  viewing  plane  which  depends  upon  the  spacecraft's 
trajectory.  The  implementation  of  this  technique  guarantees  the  maximum  information  content  of 
optical  observations  made  from  future  deep-space  probes  including  the  “Out-of-Ediptic”  mission 
scheduled  for  launch  in  1983. 


1.  INTRODUCTION 

Measurements  of  the  brightness  (Z)  and  polarization  (P) 
of  the  zodiacal  light  represent  weighted  mean  values  of 
the  volumetric  differential  scattering  cross-section,  no-,, 
where  n  is  the  number  density  of  the  dust  (cm-3)  and  tr(  is 
the  cross-section  (cm3sr-’).  The  subscript  i  refers  to  the 
two  orthogonal  components  of  die  scattered  radiation.  If 
the  electric  vector  is  perpendicular  (parallel)  to  the  scat¬ 
tering  plane,  then  i  - 1  («).  Each  elemental  volume  of 
interplanetary  space  contributes  to  the  brightness  along 
the  line -of -sight  (LOS)  in  the  amount 

&  -  ■  ”  no,  (i-J.,|),  (1) 

where  F0  is  the  tolar  spectral  irradiance 
(ergs  cm-3  s-*  A-1)  at  the  Earth,  Rra  is  the  Earth-Sun 
distance,  and  R  it  the  heliocentric  distance  to  the  volume 
element.  The  total  contribution  along  the  LOS  is  then 
given  by 

Z-Z.+Z.-f  (C^+Wdi,  (2) 

Los 

where  dl  is  a  small  increment  of  the  LOS  centered  on  the 
elemental  volume,  and 


(Henceforth,  the  subscripts  denoting  polarization  will  be 
omitted  except  when  needed  for  darity). 

The  quantity  no  or,  equivalently,  {  is  the  most  funda¬ 
mental  description  of  the  optical  properties  of  the  interp¬ 
lanetary  dust  complex;  no  optical  measurements  of  the 
zodiacal  light  can  yield  more  information  than  the  spatial, 
scattering  angle  (9)  and  wavelength  dependence  of  (.  It  is 
this  functional  dependence  which  must  be  directly  com¬ 
pared  to  light  scattering  calculations  or  experiments!  data 
to  determine  the  size,  shape  and  index  of  refraction  of 
interplanetary  dust  In  lieu  of  direct  determinations  of  {, 
various  assumptions  have  traditionally  been  employed  to 
interpret  observations  of  Z  and  P.  The  usual  practice  is  to 
assume  that  the  nature  of  the  interplanetary  dust  does  not 
depend  on  ha  location  in  the  solar  system  (<r  is  a  function 
of  9  only)  and  that  the  spatial  distribution  of  the  dust  can 
be  represented  by  a  simple  law  (e.g.  n  *  r-*  where  r  is  the 
in-ecliptic  heliocentric  distance).  While  this  practice  gives 
qualitative  insight  into  the  nature  of  the  dust  particles, 


tbeae  assumptions,  without  verification,  can  introduce 
misconceptions  about  foe  true  spatial  and  scattering  angle 
dependence  of  no. 

Dumont  (1973)  was  the  first  to  point  out  that  deep- 
space  probes  could  provide  a  direct  determination  of  £. 
thus  removing  a  major  impediment  to  understanding  the 
nature  of  the  interplanetary  dust  complex.  He  developed, 
on  physical  grounds,  a  method  of  inverting  equations  (2) 
and  (3)  for  zodiacal  light  observations  of  ( and  made  from 
within)  the  symmetry  plane.  Although  not  explicitly  stated 
in  his  paper,  his  method  relies  only  on  the  existence  of  a 
plane  of  symmetry  for  the  zodiacal  light.  Maintaining  this 
assumption  here,  we  treat  the  general  inversion  of  the 
brightness  integral  which  can  be  applied  to  an  observer 
arbitrarily  located  within  the  solar  system  and  is  particu¬ 
larly  appropriate  for  the  “Out-of-Ediptic”  mission 
scheduled  for  launch  in  1983  (Fisk  and  Axford,  1976).  In 
the  next  section,  the  nomenclature  and  geometry  neces¬ 
sary  to  evaluate  the  integrals  in  equations  (2)  and  (3)  are 
presented.  This  geometry  is  used  in  Section  3  to  derive 
foe  general  mathematical  inversion  of  those  equations.  In 
Section  4,  the  general  inversion  is  applied  to  a  single 
observer  arbitrarily  moving  through  the  solar  system.  The 
existence  of  a  unique  viewing  plane  is  established  for 
which  the  inversions  of  equations  (2)  and  (3)  are  possible. 
The  functional  dependence  of  {  on  scattering  angle  can 
then  be  determined  for  a  range  of  9  which  depends  on  the 
spacecraft  trajectory. 

2.  APPROPRIATE  GEOMETRY  AND  EVALUATION 
OF  THE  BRIGHTNESS  INTEGRAL 

Without  loss  of  generality,  the  ediptic  is  chosen  as  the 
plane  of  symmetry.  The  position  of  the  observer  (the 
point  0  in  Fig.  1)  is  defined  by  two  coordinates,  the 
ediptic  projection  (r„)  of  the  observer’s  heliocentric  dis¬ 
tance  (Rg)  and  the  observer’s  height  (hg)  above  (positive) 
or  below  (negative)  the  ediptic  plane.  Thus.  Rg3- 
r03-bh03,  and  the  angle  <ft  is  defined  by  sin  db**ha/Ro-  At 
the  point  0,  the  direction  of  the  LOS  is  defined  by  its 
ediptic  longitude  with  respect  to  the  Sun  (Ag)  and  its 
ediptic  latitude  (0).  (The  fact  that  the  plane  0  -  0  at  the 
point  0  does  not  contain  the  Sun  should  present  no 
confusion.  The  ediptic  coordinate  system  is  defined  on 
the  celestial  sphere  of  arbitrarily  large  radius,  and  thus  all 
planes  parade!  to  the  ediptic  have  0-0  whether  they 
contain  the  sun  or  not)  The  brightness  measured  at  0  can 
then  be  expressed  as  a  function  of  these  four  independent 
coordinates,  Z  -  Z(rg,  hg,  Ag,  0). 
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Flo.  2.  The  relationship  of  ths  various  angles  in  Flo. 

1  AS  DRAWN  ON  A  SPHERE  CENTERED  ON  THE  OBSERVER. 

Central  to  this  discussion  is  the  expression  (or  the 
elongation  (e),  the  Stm-observer-LOS  angle,  in  terms  of 
the  four  independent  coordinates.  Figure  2  shows  the 
relationship  of  the  various  angles  as  viewed  on  a  sphere 
centered  on  the  observer.  The  solution  of  the  spherical 
triangle,  combined  with  the  definition  of  dfe,  yields 

cos  *  “  — ^sin  0  +^“"OOs  0  cos  A£.  (4) 

Rq 

For  completeness,  we  also  list  the  expression  for  the 
indination  (i)  of  die  scattering  plane  with  respect  to  the 
ediptic, 

cos  i  =  cos  0  »  ^  cos  3.  (S) 

RgUnc  B 

where  b  and  B  are  defined  below  (Equations  8  and  9). 

The  brightness  seen  by  an  observer  at  0  is  due  to  the 
contribution  of  all  volume  elements  along  the  LOS.  An 
arbitrary  volume  element,  as  pictured  in  Fig.  1,  is  de¬ 
scribed  by  its  position  (r,  h)  and  direction  of  the  LOS 
(A',0)  in  the  same  manner  as  was  used  to  describe  the 
position  and  look  direction  of  the  observer,  but  without 
subscripts.  The  scattering  angle  (analogous  to  the  ob¬ 
server’s  elongation  angle)  is  defined  by  the  LOS  and  a  line 
(R)  from  the  sun  to  the  typical  volume  element.  The 
plane  containing  these  two  lines  is  the  scattering  plane. 
Using  relationships  similar  to  those  used  for  c,  it  can  be 
seen  that 

cos 8  =--sin  0  -t—oos  0  cos  A'  (6) 

K  K 

and  that,  for  the  volume  element  centered  on  the  ob¬ 
server,  0  =  «. 

There  exist  three  constants  which  define  any  given 
LOS.  The  angle  0  is  constant  along  any  such  line  as  are 
the  two  lengths  designated  B  and  b  in  Fig.  1.  B  is  the 
closest  solar  approach  of  the  line  of  sight,  and  b  is  the 
same  for  the  projection  of  the  LOS  in  the  ediptic.  (The 
plane  defined  by  B  and  b  is  normal  to  the  LOS  at  the 
point  Q.)  We  make  extensive  use  of  the  fact  that,  along 
any  given  LOS, 

0-  constant,  (7) 

b  "r„  sin  A£  =  r  sin  A’ »  constant,  (8) 

B  “Ro sine  -  R  sin# “ constant.  (9) 

Any  point  along  the  LOS  can  be  parameterized  by  one 


other  variable.  A  convenient  choice  is  the  acatsering 
angle,  6.  Equations  (6M9)  allow  r,  h,  and  A*  so  be 
expressed  as  functions  of  •,  0,  b,  and  B,  the  latter  three 
being  constant  along  any  given  LOS.  Explicitly, 

B  sin  0  - 

h= - -*s/Ba-»*eos0.  (10) 

tan  0 

(B  cos0  - 

-^j~±jB2-b2tmp)  ,  (11) 

b  B  cos  0  - 

- -- - — ±VB*-Vsin0.  (12) 

tan  A’  tan  6 


The  geometric  relationships  expressed  by  equations  (10)- 
(12)  may  also  be  obtained  by  considering  Fig.  1,  where 
(B’-b1)1”  is  designated  as  q.  The  +(-)  sign  in  the  * 
term  of  equation  (10)  is  used  if  die  point  Q  is  in  the 
northern  (southern)  ediptic  hemisphere.  In  equations  (11) 
and  (12),  the  positive  sign  is  used  if  0 >0  and  O  is  in  the 
northern  hemisphere  or  if  0  <0  and  Q  is  in  the  southern 
hemisphere;  the  negative  sign  is  used  otherwise. 

In  evaluating  the  brightness  (2)  and  polarization  (3) 
integrals,  it  is  convenient  to  use  the  scattering  angle  as  die 
variable  of  integration.  The  relations 

1= - —  for  »/2<e=Ssr,  (13) 

tan(*~0) 

<=—  for  0s»Ss/2  (14) 

tan  $ 


are  used  to  transform  the  LOS  integration. 


L-wz)  rj+L-.i  r’-bT* 


(15) 


(0  is  a  more  convenient  variable  of  integration  than  is  f 
because  the  latter  has  a  branch  point  at  l  =  0,  the  point  Q 
in  Fig.  1).  The  brightness  integral  then  becomes 


Z(r0,ho,  Ai,  0)“|  jV+li5){(r,M)d0.  (16) 

The  following  recipe  can  now  be  used  to  integrate  (16). 
The  position  (r0,  h0)  of  the  observer  and  his  choice  of  the 
LOS  (Ai,0)  are  used  to  determine  0,  b,  and  B  via 
equations  (7H9).  The  position  of  each  volume  element 
(r,  h)  along  the  line -of -sight  can  then  be  specified  at  a 
function  of  0  only  (0,  b,  and  B  are  constant  once  the  LOS 
is  specified)  via  equations  (10)  and  (1 1).  The  integrand  in 
equation  (16)  is  then  a  function  of  0  alone,  and  the 
evaluation  of  the  integral  it  straightforward. 


3.  INVERSION  or  THE  BRIGHTNESS  INTEGRAL 

The  inversion  of  the  brightness  integral  is  a  mathemati¬ 
cal  description  of  the  simple  physical  concept  pictured 
schematically  in  Fig.  3.  There,  die  cumbersome  three- 
dimensional  geometry  has  been  restricted  to  one  dimen¬ 
sion  so  that  the  sun  (O),  die  observer  (O),  and  line  of 
sight  (-*)  are  embedded  in  a  column  of  interplanetary 
material  of  length  x  ■  L  The  observer,  at  x  looks  in 
the  antisun  direction  <€“#“»)  and  records  a  brightness 
Z(n),  at  indicated  in  die  figure.  The  observer  then  moves 
a  small  amount  (da)  and  makes  another  measurement. 
Z(a+da).  The  difference  in  these  measurements,  dZ(«), 
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2(o)  Z  (o+do)  dZ(o) 


Flo.  3.  A  SCHEMA  TIC  REPRESENTATION  OF  THE  PHYSICAL  MEANING  OF  INVERTING  THE  BRIGHTNESS  IN¬ 
TEGRAL. 


yields  iafonnation  about  the  material  only  in  the  distance 
interval  do,  the  immediate  vicinity  o /  the  moving  ob¬ 
server.  Furthermore,  this  information  depends,  not  on  Z, 
bat  upon  its  derivative,  dZ/da. 

We  now  formally  apply  this  technique  to  equation  (16). 
The  act  of  moving  the  observer  along  the  LOS  corres¬ 
ponds  to  differentiating  the  right  hand  side  with  respect  to 
c  According  to  the  fundamental  theorem  of  calculus, 

i£[{V+h*,  K  #)<m]  -  ~  Ur*  h*  e  -  *). 

(17) 

The  left  hand  side  of  (16)  becomes. 


d  w  4Zdr,  iZ  dhg  3Z  dAg  dZd ft 

08) 

The  total  derivatives  are  geometric  quantities  which  can 
be  evaluated  from  equations  (7),  (10),  (11)  and  (12)  with 
k.  r.  A'  and  •  set  equal  to  hg.  r0,  AJ  and  «,  respectively. 
After  simplifying  the  expressions  through  the  use  of  equa¬ 
tions  (8)  and  (9),  we  And  that 


B 


cos  8  cos  A* 


(19) 


ghg  Rg*  . 

T3”  „  *“>0. 

da  B 

(20) 

^-^cosfisinA;. 

(21) 

*-a 

d« 

(22) 

The  inversion  of  the  brightness  integral  is  therefore  given 
by 

B  dZ  I 

«r*ho.«-0)--R#>— ^ 


chosen  to  be  in  the  direction  or  antidirection  of  motion, 
two  measurements  of  Z  separated  by  an  interval  de  along 
the  LOS  immediately  provide  Ur*  K,,*m  9).  Thus,  C  can 
always  be  determined  for  two  values  of  9  (differing  by  w) 
at  all  points  along  the  spacecraft  trajectory.  The  func¬ 
tional  dependence  of  {  on  9(r,  h.  A’,  ft)  is  contained  in  the 
second  statement  of  equation  (23),  but  it  is  not  possible 
for  a  tingle  spacecraft  to  measure  all  three  partial  deriva¬ 
tives  for  even  a  tingle  point  in  the  solar  system.  A  moving 
observer  can,  at  most,  measure  the  angular  derivative  and 
one  spatial  derivative  (in  the  direction  of  motion).  Dum¬ 
ont  et  aL  (in  presl)  have  demonstrated  that  those  two  meas¬ 
urements  are  sufficient  to  invert  the  brightness  integral  for 
directions  confined  to  a  special  viewing  plane  referred  to 
here  as  the  plane  of  inversion.  This  plane  limits  the  range 
for  which  the  functional  dependence  of  (('.  K  9)  on  9  can 
be  derived  from  observations.  In  the  mathematical  de¬ 
velopment  which  follows,  the  uniqueness  of  the  viewing 
plane  is  established  and  the  inversion  formula  it  cast  into 
a  form  similar  to  that  used  for  in -ecliptic  observations. 

Let  the  unit  vector  l  represent  a  segment  of  the  space¬ 
craft  trajectory  in  the  direction  of  motion;  i.e.  the  space¬ 
craft's  instantaneous  speed  is  ds/d(.  The  plane  that  contains 
i  and  the  normal  to  the  r,  h  plane  is  the  plane  of  inver¬ 
sion.  Let  4  be  the  angle  between  the  inversion  plane  and 
the  ecliptic  and  define  a  new  coordinate  system  by 

x~rootd  +  hsln<5, 

y--rsind  +  hoos4,  (24) 

as  shown  in  Fig.  4.  The  two  spatial  derivatives  in  equation 
(23)  then  transform  via 

4Z_4Ztas+8Z»yg 
llg  9Xq  if q  iy0  hp 

^  4Z  d*p+£Z  «5) 

lip  Itplip  dypdhg 


so  that 


9Z 


az 


az 


»  —  coaffeoaAi-  —  sin? - —  eos0sinA;.  (23) 


*8o 


r.*A. 


az 

Ur*  *o.  •  “  «) "  —  (coa  ft  cos  A  cos  d  -sin  ft  sin  4) 

**o 


The  polarization  per  unit  volume  ■  obtained  by  applying 
equation  (23)  to  each  component  of  Z,  separately  and 
then  forming  the  quantity  For  the  case  ft  -0, 

h,“0,  the  general  solution  (23)  reduces  to  the  two- 
dimensional  (in-ecliptic)  form  given  by  Leinert  (1975V 


az 

- (cot  ft  cos  A i  tin  d + sin  ft  cos  d) 

*7e 


az 


(26) 


T-cospsinAg. 


The  coefficient  of  aZfay0  in  equation  (26)  vanishes  for 
4.  APPLICATION  TO  OUT-Or-ECUTOC  ‘hose  viewing  directions  in  which 


OBanVATTONB 

The  implementation  of  equation  (23)  requires  the  ob- 
aerver(s)  to  be  moving  through  the  solar  system.  The  fttst 
statement  of  the  equation  formalizes  the  descriptive  argu¬ 
ment  given  at  the  beginning  of  Section  3.  If  the  LOS  is 


oosAi- 


rang 

rand 


(27) 


Equation  (27)  is  the  coordinate  definition  of  the  plane  of 
inoenkm. 


i 
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Flo.  4.  Construction  of  the  inversion  (shaded)  plane-,  it  is  normal  to  the  r,  h  plane  and  con¬ 
tains  S,  THE  DIRECTION  OF  MOTION  OF  THE  SPACECRAFT. 

The  angle  <t>  is  the  inclination  of  the  inversion  plane  with  respect  to  the  ecliptic. 


For  the  case  of  <5  =  0,  it  is  obvious  from  the  coefficient  of 
3ZJdy0  in  equation  (26)  that  the  inversion  plane  is  parallel 
to  the  ecliptic;  i.e.  0  =  0. 

For  viewing  directions  confined  to  the  inversion  plane,  £ 
can  be  determined  without  knowledge  of  3ZJ3 y„,  pro¬ 
vided  that  3ZJ3xa  and  3ZJ3 A„  are  known.  The  angular 
derivative,  3ZJ3\'0,  is  obtained  in  practice  by  scanning 
across  the  plane  of  inversion  in  directions  parallel  to  the 
ecliptic.  The  spatial  derivative,  3ZJ3xa,  is  obtained  by 
differencing  two  sets  of  observations  which  are  separated 
by  a  small  distance  in  x,  where  8x0  is  the  projection 
(x  =s  cos  v  where  v  is  measured  in  the  inversion  plane  as 
in  Fig.  4)  of  8s  on  the  x  axis.  It  is  important  to  recognize 
that  the  spacecraft  need  not  be  traveling  parallel  to  the  x 
axis  in  order  to  determine  3ZJ3x0.  This  partial  derivative 
is  a  true  gradient:  the  maximum  rate  of  change  of  Z  is 
along  the  x  axis.  Because  of  the  assumption  of  cylindrical 
symmetry,  this  gradient  has  a  zero  component  in  the 
direction  normal  to  the  r,  h  plane,  and  its  value  is  im¬ 
mediately  determined  from  increments  in  Z  and  x  only  or 


Flo.  5.  The  obsf.rver-s  viewing  sphere  showino  the 
INVERSION  (SHADED)  PLANE. 


by  noting  that 


8Z  3Z 

—  =  —  COS  V. 

3s  3x0 


(28) 


The  non-vanishing  coefficients  in  equation  (26)  have  an 
obvious  geometric  interpretation  as  shown  on  the  ob¬ 
server’s  viewing  sphere  in  Fig.  5.  If  the  angle  £  is  meas¬ 
ured  in  the  inveision  (shaded)  plane  from  the  solar  meri¬ 
dian  to  the  viewing  direction,  then  the  solution  of  the 
appropriate  spherical  triangle  in  Fig.  5  yields 

cos  £  =  cos  0  cos  A^cos  <(>  -sin  0  sin  <f>,  (29) 

sinf  =  cos0  sin  (30) 

and  t  and  £  are  related  by 

cose  =  cos  £  cos  (ij>0-i5).  (31) 

These  substitutions  bring  the  inversion  formula  into  the 
simple  form 

£(r0,h  o-0  =  «)  =  7— «»£ — “~r  sin  (•  (32) 

3x0  r08Ao 

which  is  reminiscent  of  the  in-ecliptic  solution  and  re¬ 
duces  to  it  (x0=r0,  £  =  A„  =  t)  for  $  =  0,  0  =  0,  ^,=0. 
Formula  (32)  allows  £  to  be  determined  at  all  points  of 
the  solar  system  encountered  by  the  probe  and  for  all 
values  of  the  scattering  (elongation)  angle  within  the 
range 

(<jr0-(5)se  =  €Sw-(i|ro-<5).  (33) 

The  determination  of  £  requires  observations  to  be 
made  in  directions  other  than  the  inversion  plane.  The 
angular  derivative  in  equation  (32)  can  only  be  obtained 
from  observations  which  at  least  straddle  this  plane.  The 
fact  that  8Z/8 A0’,  rather  than  8Z/8£,  appears  in  equation 
(32)  is  due  once  again  to  the  only  crucial  assumption 
inherent  in  this  development,  that  of  cylindrical  sym¬ 
metry. 

A  very  direct  check  on  this  sole  assumption  wffl  be 
possible  during  the  Out-of-Ecliptic  mission.  After  a  gravi¬ 
tational  assist  by  Jupiter,  the  spacecraft  will  travel  back 


SS6 


Reaearch  note 


“over”  (“under”)  the  Sun.  Approximately  3  yr  into  the 
minion,  solar  polar  passage  wOI  be  attained  at  which  time 
the  spacecraft  will  be  near  the  center  of  the  solar  system 
but  approximately  1.6  a.u.  above  (below)  the  ediptic 
plane.  From  tins  vantage  point,  all  isophotes  of  the  bright¬ 
ness  of  the  zodiacal  light  will  be  circles  centered  on  the 
sun  if  the  assumption  of  cylindrical  symmetry  is  correct. 
Any  long-term  influences  of  the  planets  on  the  interp¬ 
lanetary  dust  complex  should  also  be  discernaUe  as 
changes  in  the  brightness  distribtuion  near  the  orbits  of 
the  planets. 
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increase  in  the  visual  luminosity,  with  little  change  in  the 
bolometric  luminosity.  Outflow  velocities  will  be  lower  than  in 
classical  novae:  roughly  the  same  order  as  the  escape  velocity 
from  the  stellar  surface. 

Consequently  there  will  be  two  types  of  spectroscopic  and 
colour  change  produced  by  variations  in  the  accretion  rate.  In 
the  subcritical  regime,  increased  accretion  rates  produce 
increasing  disk  temperatures  along  with  increasing  luminosity — 
both  visual  and  bolometric.  In  the  supercritical  regime  the  disk  is 
surrounded  by  an  outflowing  optically  thick  wind  and  the 
bolometric  juminosity  is  held  close  to  the  Eddington  value, 
Increases  in  the  wind  mass-loss  rate  due  to  further  increased 
accretion  now  produce  decreasing  temperatures  along  with 
increasing  visual  luminosity.  The  emission  line  spectrum  of  the 
disk  will  be  replaced  by  an  absorption  line  spectrum  formed  in 
the  outflowing  envelope — similar  to  the  F  or  A  type  spectrum 
observed  in  classical  novae  at  outburst  maximum. 

Is  there  any  evidence  for  this  supercritical  state?  The  photo¬ 
graphic  magnitude  of  M33  Var  A  fell  by  3.S  mag  (ref.  5)  in  the 
visual  in  1950-52,  rose  by  1.0  mag  in  1952,  before  finally  falling 
off  sharply.  Hubble  and  Sandage  note  that  the  changes  in  colour 
index  indicate  that  the  bolometric  luminosity  changed  by  only 
1.4  mag.  The  visual  luminosity  fell  by  —60  while  the  bolometric 
luminosity  fell  by  —4.  This  agrees  with  the  expected  changes  due 
to  a  variable  optically  thick  wind  at  Eddington  luminosities. 

t)  Carina  is  a  close  galactic  relative  to  the  Hubble-Sandage 
variables.  It  has  broad  spectral  lines,  with  relatively  strong  Fc  n 
and  [Fe  n]4.  Davidson*9  deduced  a  temperature  of  3  x  ICr  K  for 
the  underlying  star.  Between  1836  and  1856  17  Carina  under¬ 
went  a  major  outburst,  followed  by  a  subsidiary  outburst  in 
1889.  The  early  F  absorption  spectrum  of  1893  was  slowly 
replaced  by  the  characteristic  emission  spectrum.  The  absolute 

magnitude  at  outburst  was - 14  (ref.  20),  several  magnitudes 

brighter  than  normal  classical  novae.  At  present,  a  cool  sur¬ 
rounding  dust  cloud1'  heated  by  the  underlying  stellar  source 
produces  the  bulk  of  the  emission.  The  bolometric  luminosity  is 
still  —  -14,  but  most  of  the  emission  is  re-radiated  IR  from  the 
dust  shell.  This  high  and  constant  luminosity,  the  spectral 
changes,  the  broader  than  average  emission  lines,  and  the 
nova-like  explosion  indicate  that  17  Carina  could  contain  a 
massive,  main-sequence,  accreting  component  in  a  binary  which 
underwent  a  burst  of  supercritical  accretion.  Polarisation 
observations10  of  the  nebula  suggest  that  the  dust  is  concen¬ 
trated  in  an  annulus  of  size  104  au,  and  this  has  been  interpreted 
as  evidence  of  binary  structure.  This  is  the  first  observational 
indication  of  possible  binary  structure  in  this  type  of  object.  A 
search  for  binary  eclipses  in  Hubble-Sandage  variables  with 
periods  in  the  range  weeks  to  years  would  test  the  binary  model. 

The  bright  galactic  source  MWC349  is  a  similar  object7. 
Thompson  et  al l.11  interpret  optical  and  IR  observations  as  a 
proto-planetary  disk  seen  face-on.  In  terms  of  the  structure  of 
the  luminous  source  their  model  closely  resembles  that 
described  here.  The  main  difference  is  their  interpretation  as  a 
preplanetary  disk  rather  than  an  interacting  binary. 

If  the  brightest  known  optical  stars  are  powered  by  infall  of 
material  on  to  a  main-sequence  star,  the  importance  of  accretion 
as  an  astrophysical  energy  source  would  be  established  beyond 
its  present  applications  to  X-ray  binaries,  novae  and  active 
galactic  nuclei  and  quasars11.  Such  binaries  should  exist,  and 
they  would  be  classified  as  Hubble-Sandage  variables  when 
observed  in  external  galaxies. 

I  thank  J.  Barrow  and  A.  Ed'/ards  for  helpful  comments. 
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Deep-space  probes  offer  excellent  opportunities  to  study  inter¬ 
stellar  grains  streaming  into  the  Solar  System.  Direct  detection 
of  such  particle*  would  considerably  add  to  our  knowledge  in 
many  areas  of  cosmic  dost  research  and  could  also  provide 
information  on  the  characteristics  of  both  the  interstellar 
medium  and  the  interacting  solar  wind.  Here  we  present  the 
results  of  a  theoretical  study  of  the  interaction  between  inter¬ 
stellar  grains  and  the  solar  wind,  including  the  effects  of  solar 
cycle  variations  in  the  interplanetary  magnetic  field.  The  latter  is 
shown  to  influence  significantly  regions  of  concentration  and 
depletion  of  interstellar  grains  within  the  Solar  System. 

Observations  of  the  neutral  gas  in  the  solar  vicinity  suggest 
that  the  Solar  System  is  moving  relative  to  the  local  interstellar 
matter  in  a  direction  =7°  north  of  the  solar  equatorial  plane 
(ecliptic  longitude  and  latitude:  A  =*250°.  /3“7°)*,  and  at  a 
speed  of  “20  km  s~*  (ref.  2).  Bertaux  and  Blamont1  computed 
number  densities  of  interstellar  grains  in  the  solar  cavity  taking 
into  account  only  gravitation  and  radiation  pressure.  Of  the 
types  of  grains  used  in  their  model  and  based  on  their  calculation 
of  p  (radiation  pressure  to  gravitational  attraction)  only  0.02- 
lim  radius  silicate  grains  were  able  to  penetrate  deeply  inside  the 
Solar  System.  They  estimated  the  number  density  of  these  dust 
grains  to  be  strongly  peaked  on  a  line  extending  from  the  Sun 
backwards  in  the  downstream  direction.  This  number  density 
was  found  to  be  higher  by  a  factor  of  —100  than  direct  obser¬ 
vations  of  dust  grains  by  Alvarez4.  In  an  attempt  to  clarify  this 
difference.  Levy  and  Jokipii5  made  an  order  of  magnitude 
estimate  of  the  Lorentz  force  exerted  by  the  solar  wind  magnetic 
field  on  a  charged  grain.  They  showed  that  if  the  0.02-pm  radius 
silicate  grains  of  Bertaux  and  Blamont  carry  a  net  charge  of  5  V, 
then  the  Lorentz  force  strongly  dominates,  thereby  destroying 
the  gravitational  focusing  and  suggesting  that  these  grains  are 
largely  excluded  from  the  Solar  System. 

Studies  of  interstellar  grains  suggest  two  populations  of  parti¬ 
cles*:  (1)  elongated  core-mantle  particles  of  characteristic 
radius  “0.12  pm,  with  a  “0.07  pm  radius  core  of  silicate  type 
material,  and  a  rather  complex  mantle  of  molecular  mixtures  of 
mainly  oxygen,  carbon  and  nitrogen;  and  (2)  “0.01  pm  thick 
silicate  and/or  graphite  type  grains.  Bertaux  and  Blamont  did 
not  consider  any  of  these  populations.  In  the  following  cal¬ 
culations  (see  also,  ref.  1 3)  we  show  that  although  population  ( 1 ) 
particles  have  p  >  1,  they  still  penetrate  deep  inside  the  Solar 
System. 

The  magnitude  of  the  radiation  pressure  is  given  by  Pp,  = 
tS/c )0p,,  where  S  is  the  flux,  c  is  the  speed  of  light,  and  Q.,  is  the 
radiation  pressure  efficiency:  Q„  =  G.«,+  C>„.  (1  —  (cos  <?)). 
where  0<hI  and  Q„.  are  the  efficiency  factors  for  absorption  and 
scattering,  respectively,  and  (cos  0)  is  the  asymmetry  factor  over 
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Fig.  1  Instantaneous 
positions  of  the.  0.12  jim 
core-mantle  particles 
initially  moving  in  the 
direction  A  ■=  70’,  0  =  7“ 
along  the  solid  lines.  The 
particles  are  entering  the 
Solar  System  at  X  = 

+50  au  with  2  yr  interval. 

Three  streams  (■,  x,  •) 
introduced  at  Z  =  +6  au 
are  projected  onto  the 
solar  equatprial  (X -Y) 
plane  at  the  beginning  (a) 
and  in  the  middle  (cTof  a  T 
solar  magnetic  cycle.  In  b  ■  ■ 
and  d  the  three  streams  — * 
(0,  x,  ©)  introduced  at 
y *=0  are  similarly  pro¬ 
jected  oglo  a  plane  (X-Z) 
perpendicular  to  the  solar 
equatorial  plane  (dotted 
line). 


the  scattering  angle  0.  Scattering  dominates  over  absorption  for 
these  dielectric  particles  in  the  scattering-efficient  resonance 
region,  and  Op,  was  estimated  using  Mie  theory  for  homo¬ 
geneous  spheres  of  equal  volume  and  averaged  core-mantle 
index  of  refraction.  For  a  spherical  silicate  core  coated  by  a 
concentric  smooth  spherical  ice  mantle  representing  particles  of 
population  (1)  the  computed  Op,  value  using  this  approximation 
was  found  to  be  0.4  (0  =  1.3).  As  the  grains  that  penetrate 
(population  (1))  spend  30  yr  or  less  in  the  Solar  System,  all 
radiation  forces  other  than  radiation  pressure  arc  negligible. 
Both  solar-induced  photoelectron  emission  and  the  accretion  of 
solar  wind  ions  will  tend  to  charge  these  grains  positively, 
whereas  electron  accretion  produces  negative  charge.  The 
equilibrium  charge  state  is  given  by  the  relation 

<r.(V0)f.-<rl(Va)fi  +  ycr,hf,h  (1) 


where  the  tr s  are  the  cross-sections  of  the  grains  for  U  V  photons, 
ions  and  electrons,  the  /s  are  the  corresponding  fluxes,  y  is  the 
photoelectron  yield,  and  V0  is  the  electrical  potential  of  the 
grains.  <r  and  y  have  been  estimated  for  numerous  silicates,  and 
these  lead  to  charges  ranging  from  5  to  12  V  (refs  7,  8).  The 
photoemission  properties  are  not  usually  decreased  by  surface 
impurities  or  adsorbed  layers7,  or  by  amorphous  structure*. 
However,  for  ice  the  work  function  is  near  the  Lyman  edge,  and 
therefore  photoemission  is  less  efficient.  We  will  adopt  here 
+7  V  for  bare  particles,  and  +1.2  V  for  the  core-mantle  grains. 
As  ihe  charged  grains  interact  with  the  solar  wind  described  by 
Svalgaard  el  a/.*,  they  will  be  subject  to  a  Lorentz  force  given  by 

frL  =  J[(V,xB)-(V,.xB)]  (2) 


where  q  is  the  grain's  charge,  V,  and  V„  are  the  speeds  of  the 
grain  and  of  the  solar  wind,  respectively,  c  is  the  speed  of  light 
and  B  is  the  magnetic  field  strength. 

The  magnetic  held  is  torn  in  an  archimedian  spiral  by  the  solar 
rotation.  Thus 


B, 


B9 

■j 


and  B*  = 


-B0wfr!(r-fi)  sin  8 

vZT1 


(3) 


where  B0  is  the  magnetic  held  strength  at  some  reference 
distance  ( b  =  5-10  solar  radii),  o>  is  the  solar  angular  velocity,  r  is 
the  heliocentric  distance,  <t>  is  the  azimuth  angle,  and  0  is  the 


heliocentric  colatitude  angle.  The  second  (electrical)  part  of 
equation  (2)  which  acts  along  0  will  dominate  as  the  grains  enter 
the  Solar  System  at  low  and  medium  heliocentric  latitudes.  This 
force,  acting  on  a  grain  of  given  charge,  is  independent  of  the 
speed  of  a  radially  streaming  solar  wind  as  can  be  seen  by 
substituting  equation  (3)  in  equation  (2);  high  velocity  plasma 
streams  will  not  alter  our  results  considerably. 

The  positive  and  negative  magnetic  fields  are  separated  by  a 
sinusoidal  warped  neutral  sheet10-",  the  amplitude  of  which  we 
assume  to  decay  exponentially  with  increasing  heliocentric  dis¬ 
tance.  The  speed  of  the  solar  wind  as  a  function  of  8  was  assumed 
to  be  of  the  form 


V„~a  +  b  |cos  0|  (4) 

where  a  is  the  mean  solar  wind  velocity  at  the  equator 
(400  km  s'*)  and  b  is  300  km  s~\ 

The  grains  also  encounter  other  less  important  forces  such  as 
Coulomb  drag  due  to  the  nearby  passing  solar  wind  electrons, 
which  is  given  by 


vs'N.e'Vl 

m.W1 


In 


(5) 


where  s  and  V0  are  the  radius  and  charge  of  the  grain,  respec¬ 
tively,  N„  e,  and  m,  are  the  electron  number  density,  charge  and 
mass,  respectively,  and  W  is  the  velocity  of  the  grains  relative  to 
the  electrons.  P  is  the  impact  parameter,  which  is  assumed  to  be 
“ Nl /J,  following  a  poissonian  distribution  for  the  nearest 
neighbour  distance.  Another  force  is  ion  pressure  resulting  from 
impinging  solar  wind  ions  on  the  grains  which  is  given  by 

F,  =  WM,VF  (6) 

where  F  is  the  flux  of  the  solar  wind  ions,  V  is  their  velocity  and 
Mi  is  the  mean  ion  mass  =2x1 074  g. 

Trajectories  of  interstellar  grains  streaming  into  the  Solar 
System  were  computed  as  follows:  the  grains  are  introduced  at 
the  boundary  of  the  solar  wind  cavity  (X  =  50  AU)  at  the  begin¬ 
ning  of  the  22-yr  solar  magnetic  cycle  (defined  as  the  cycle 
beginning  when  the  solar  dipole  field  acquires  a  south 
component).  The  equations  of  motion  are  integrated  at  small 
intervals  of  time  (24  h),  and  the  position,  speed  and  acceleration 
of  the  particles  are  determined. 


Distance  (au)  X 

Fig.  2  Instantaneous  positions  of  combined  high  latitude  streams  and  the  streams  shown  in  Fig.  lb  and  d,  at  time  zero  (a )  and  time  1 1  yr  (6)  of  a 

solar  magnetic  cycle. 


The  results  of  these  computations  are  shown  in  Figs  1  and  2. 
Figure  la  and  c  shows  the  projection  onto  the  solar  equatorial 
plane  (XT-plane)  of  the  instantaneous  positions  of  population 
particles  arising  from  three  streams  which  differ  only  by  a 
separation  of  10  AU  in  the  initial  value  of  their  Y  coordinate. 
Figure  lb  and  d  shows  the  instantaneous  positions  in  the 
XZ-plane  of  particles  arising  from  three  streams  with  Y= 0. 
One  stream  is  initially  directed  towards  the  Sun  and  the  other 
two  lie  parallel  to  and  10  AU  on  either  side  of  that  direction.  The 
sequence  of  points  belonging  to  a  stream  corresponds  to  posi¬ 
tions  of  particles  which  differ  by  2  yr  from  the  initial  time  when 
each  particle  was  introduced  at  the  boundary  of  the  solar  wind 
cavity.  Figure  la  and  b  shows  the  instantaneous  positions  of  all 
the  particles  at  a  time  corresponding  to  the  beginning  of  the  solar 
magnetic  cycle.  Figure  lc  and  d  repeat  this  information  for  a 
time  1 1  yr  into  the  solar  magnetic  cycle.  The  contrast  between  a, 
b  and  c,  d  clearly  demonstrates  the  decisive  influence  of  the  solar 
magnetic  cycle  on  the  distribution  of  the  grains.  The  streams  are 
diverted  towards  the  X-axis  in  Fig.  la  and  away  from  that  axis  in 
Fig.  lc.  Figure  16  shows  that  when  the  solar  dipole  field  axis 
points  north,  the  grains  are  diverted  from  the  solar  equatorial 
plane,  mainly  due  to  the  action  of  the  electrical  part  of  the 
Lorentz  force.  The  grains  are  concentrated  towards  this  plane 
when  the  solar  dipole  field  points  south  (Fig.  Id).  Some  3  to 
4  AU  ahead  of  the  Sun  the  streams  cross  the  equatorial  plane, 
creating  an  enhancement  in  the  number  density.  The  grains 
behind  the  Sun  experienced  the  oppositely  directed  magnetic 
field  of  the  first  half  of  the  solar  magnetic  cycle  as  they  entered 
the  solar  wind  cavity,  thus  creating  a  depletion  zone  in  the  wake 
of  the  Sun. 

Our  computations  are  almost  insensitive  to  whether  the 
reversal  of  the  direction  of  the  interplanetary  magnetic  field  is 
instantaneous  or  whether  it  is  as  gradual  as  a  sine  function. 
These  results  are  based  on  a  simple  dipole  field  where  higher 
order  multipole  or  asymmetrical  magnetic  fields  are  ignored. 

In  Fig.  2a  and  b  additional  streams,  projected  in  the  XZ 
plane,  are  introduced  higher  in  latitude  than  the  streams  in  Fig. 
16  and  d.  Figure  2a  shows  a  concentration  of  particles  develo¬ 
ping  along  a  line  X  =  -20  AU,  Z  =  ±  10-20  AU  (peaked  around 
<b  =  180°).  There  is  also  a  depiction  zone  of  radius  =20  AU 
surrounding  the  Sun.  The  deviation  from  the  XZ-plane  is  =£  ± 
0.1  AU.  The  value  of  Qp(  used  in  the  calculation  is  0.4,  which 
might  be  increased  by  surface  roughness;  however,  this  value 
decreases  for  elongated  particles.  The  overall  effect  of  increas¬ 
ing  the  radiation  pressure  efficiency  Qp,  of  the  grains  from  0.4  to 
as  much  as  1 .0  (/)  =  3.25)  is  to  move  the  concentration  in  Fig.  26 


=3  AU  along  the  positive  X-axis.  On  the  other  hand,  a  decrease 
in  the  value  of  (?„,  has  a  negligible  effect.  The  smaller  bare  grains 
that  enter  at  low  and  medium  heliocentric  latitudes  are  essen¬ 
tially  trapped  around  the  solar  wind  magnetic  field  lines,  and 
they  are  swept  out  near  the  solar  wind  cavity. 

We  have  extended  the  study  of  the  dynamics  of  interstellar 
grains  to  particles  having  radii  of  0.12  and  0.005  p.m  streaming 
into  the  Solar  System  and  included  the  effects  of  the  solar  dipole 
field's  change  of  direction  during  the  sole  magnetic  cycle. 
Additional  forces  on  the  particles  such  as  Coulomb  drag  and 
solar  wind  ion  pressure  are  less  important  than  the  Lo'cntz 
force.  Particle  trajectory  computations  show  that  an  enhance¬ 
ment  of  the  order  of  50  or  more  in  the  number  density  of 
interstellar  core-mantle  particles  may  occur  some  3-4  au  ahead 
of  the  Sun  and  close  to  the  solar  equatorial  plane  during  yr  8-17 
of  the  solar  magnetic  cycle.  During  the  rest  of  the  cycle,  concen¬ 
trations  develop  some  20  AU  away  from  the  Sun,  while  a 
depletion  zone  of  =20  au  in  radius  surrounds  the  Sun. 

As  particle  configuration  varies  strongly  with  the  solar 
magnetic  cycle,  its  effect  must  be  taken  into  account  in  any 
treatment  of  this  problem.  Note  that  even  grains  with  p  >  1  will 
penetrate  deep  inside  the  Solar  System  during  part  of  the  solar 
magnetic  cycle. 

Subsequent  computations  may  predict  more  precisely  the 
directions  and  times  to  search  for  an  interstellar  component  of 
the  zodiacal  cloud.  Note  that  during  the  Solar  Polar  Mission 
(mid  and  late  1980s),  concentrations  of  interstellar  grains  are 
expected  (Figs  Id,  26)  which  will  offer  a  unique  opportunity  to 
detect  the  presence  of  interstellar  grains  in  the  Solar  System”. 
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ABSTRACT 

Traditionally,  earth-based  observations  of  the  zodiacal  light  (ZL) 
require  two  assumptions  for  further  analysis:  (Al)  the  dust  density  (n) 
is  a  power  of  heliocentric  distance  (R),  n  R“v;  (A2)  the  nature 
(scattering  cross  section,  a)  of  the  dust  is  independent  of  location, 
o(r,h,G)=o(0) .  Observations  from  Pioneer  10  do  not  verify  these  assump¬ 
tions. 


Both  two-dimensional  inversion  techniques  (Schuerman,  1979a)  and 
more  traditional  methods  of  analysis  (Hanner  et  al.,  1976)  have  been 
performed  on  portions  of  the  ZL  observations  made  from  Pioneer  lO/ll 
spacecraft.  In  developing  the  three-dimensional  form  of  the  ZL  inversion 
(Schuerman,  1979b;  henceforth  called  Paper  I),  it  was  realized  that  this 
more  general  technique  should  be  applied  to  the  Pioneer  data  because  the 
velocity  vectors  of  the  spacecraft  were  not  confined  to  the  ecliptic 
plane.  The  previously  mentioned  analyses  examined  data  with  viewing 
directions  parallel  to  the  ecliptic;  the  three-dimensional  inversion 

requires  viewing  along  a  great 
c*rc*e  (see  Fig*  1)  containing 
Spacecraft — j  CO  the  line  (n)  which  is  both 

/  !  normal  to  the  observer-Sun 
/*  line  and  parallel  to  the 
*R  •/  ecliptic  and  (2)  the  space- 

X/  " * - >  craft  velocity  vector  (S). 

This  great  circle  is  called 

Sun  - r - A  the  inversion  plane  in  Paper 

y  /  I  where  its  geometric  con- 

/  Ecliptic  struction  and  the  notation 

- - - •  used  here  are  defined  in  more 

Fig.  1.  Definition  of  the  inversion  detail.  The  angle  $  is  the 

(shaded)  plane;  it  is  normal  to  the  inclination  of  the  inversion 

r,h  plane  and  contains  S  ,  the  direction  plane  to  the  ecliptic,  and 
of  motion  of  the  spacecraft.  viewing  in  the  direction  of  n 


— 

. y 


/  Ecliptic  • 

Fig.  1.  Definition  of  the  inversion 
(shaded)  plane;  it  is  normal  to  the 
r,h  plane  and  contains  S ,  the  direction 
of  motion  of  the  spacecraft. 
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ABSTRACT 

Analysis  of  zodiacal  light  observations  from  Mt.  Haleakala,  Hawaii 
show  that  the  symmetry  plane  of  the  zodiacal  cloud  near  1  A.U.  is  close 
to  the  invariable  plane  of  the  solar  system.  Since  the  symmetry  plane 
of  the  inner  zodiacal  cloud  is  close  to  the  orbital  plane  of  Venus 
(Misconi  and  Weinberg,  1978;  Leinett‘  et  al.>  1979),  we  suggest  that  the 
symmetry  plane  changes  inclination  with  heliocentric  distance. 


INTRODUCTION 

The  location  in  space  of  the  plane  of  maximum  dust  density  (symmetry 
plane)  of  the  zodiacal  cloud  has  important  implications  for  the  dynamics 
of  solid  particles  in  the  solar  system.  These  implications  are  related 
to  questions  such  as  the  gravitational  perturbations  of  the  planets  on 
the  dust,  the  effects  of  solar  magnetic  forces,  and  the  general  dynamical 
history  of  dust  ejected  from  comets.  Helios  1  and  2  zodiacal  light 
measurements  near  the  region  of  1  A.U.  (Leinert  et  at 1979)  find  that 
the  symmetry  plane  in  this  region  has  an  inclination  i  *»  3  ±  0.3'’  and 
ascending  node  ft  **  87  ±  4°,  which  is  different  from  that  of  the  invariable 
plane  of  the  solar  system  (i  =  1.6®,  fl  =  107®).  At  the  same  time,  ground- 
based  observations  at  Tenerife  in  the  Canary  Islands  (Dumont  and  Sanchez, 
1968)  and  satellite  observations  from  the  satellite  D2B  (Dumont  and 
Levasseur-Regourd ,  1978)  find  the  symmetry  plane  close  to  the  invariable 
plane  near  1  A.U. 

In  the  present  study  an  attempt  is  made  to  obtain  additional  data 
on  the  location  of  the  symmetry  plane  near  the  Earth's  orbit  using 
ground-based  observations  of  the  zodiacal  light  from  Mt.  Haleakala, 

Hawaii  taken  by  Weinberg  and  Mann  (1967). 


THE  OBSERVATIONS 

Two  nights  of  observations  were  found  to  be  suitable  for  this 
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ABSTRACT 

In  an  earlier  paper  Sparrow  et  at.  (1976)  found  the  polarized  bright¬ 
ness  of  zodiacal  light  to  have  solar  color  at  five  sky  positions  for  which 
there  were  fixed-position  observations  from  Skylab:  north  celestial  pole, 
south  ecliptic  pole,  vernal  equinox,  and  two  places  near  the  north  galac¬ 
tic  pole.  The  brightness  and  degree  of  polarization  of  zodiacal  light  at 
these  sky  positions  are  derived  using  Pioneer  10  observations  of  back¬ 
ground  starlight  from  beyond  the  asteroid  belt  (Weinberg  et  at.,  1974; 
Schuerman  et  at.,  1976)  and  the  assumption  that  the  zodiacal  light  is 
also  solar  color  in  total  light. 


Ten-color  observations  of  sky  brightness  and  polarization  were  made 
of  portions  of  the  antisolar  hemisphere  from  Skylab  (Weinberg  et  at., 

1975;  Sparrow  et  at.,  1977).  Three  independent  quantities  were  determined 
directly  from  the  measurements:  total  brightness  (Bt),  polarized  bright¬ 
ness  (Bp) ,  and  orientation  of  the  plane  of  polarization  (y)  >  the  last  more 
precisely  referred  to  as  azimuth  of  vibration.  These  three  quantities  are 
related  to  the  Stokes  parameters  I,  Q,  and  U,  respectively. 

Except  in  regions  relatively  close  to  the  Milky  Way,  the  observed 
direction  and  amount  of  polarization  (x,Bp)  in  these  space  measurements 
can  generally  be  attributed  entirely  to  zodiacal  light.  To  isolate  the 
zodiacal  light  brightness,  it  is  necessary  to  remove  the  brightness  con¬ 
tribution  of  all  discrete  or  ’’resolved"  stars  in  each  6“  diameter  field 
of  view  (FOV)  and  to  subtract  the  background  starlight  (integrated  star¬ 
light,  diffuse  galactic  light,  light  from  extragalactic  sources).  Thus, 
the  zodiacal  light  total  brightness  (Bz)  is  obtained  as  the  residual  after 
a  several-step  subtraction  process.  This  subtraction  process  has  been  the 
largest  single  source  of  error  in  determining  the  brightness  of  zodiacal 
light  in  observations  from  space  and,  even  more  so,  in  observations  from 
the  ground . 

Skylab  observations  included  sky-scanning  and  fixed-position  measure¬ 
ments,  the  latter  at  the  north  celestial  pole,  south  ecliptic  pole,  and 
vernal  equinox,  and  at  two  positions  near  the  north  galactic  pole.  A 
number  of  sequences  of  filters  (Table  1)  were  used  to  observe  these 
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ABSTRACT 

The  classical  restricted  three-body  problem  is  generalized  to  in¬ 
clude  the  force  of  radiation  pressure  and  the  Poynting-Robertson  effect. 
The  positions  of  the  Lagrangian  points  L^  and  L^  are  found  as  functions 
of  P>  the  ratio  of  radiational  tp  gravitational  forces.  The  Poynting- 
Robertson  effect  renders  the  L^  and  L5  points  unstable  on  a  time  scale 
(T)  long  compared  to  the  period  of  rotation  of  the  two  massive  bodies. 
For  the  solar  system,  T  is  given  by  T  *=  [(l-g)^/g)  544  a2  yr  where  a 
is  the  separation  between  the  Sun  and  the  planet  in  AU.  Implications 
for  space  colonization  and  a  mechanism  for  producing  asymmetries  in  the 
interplanetary  du&t  complex  are  discussed;  testing  the  latter  may  be 
possible  from  the  Zodiacal  Light/Background  Starlight  Experiment  aboard 
the  International  Solar  Polar  Mission  spacecraft  to  be  launched  in  1983. 


A  steady  state  solution  of  the  motion  of  a  p  particle  (a  particle 
whose  surface  to  volume  is  such  that  the  force  of  radiation  pressure  is 
non-negligible)  in  the  gravitational  field  of  two  larger,  orbiting 
bodies  of  masses  and  M2  is  sought.  It  is  assumed  that  Mj  and  M2  are 
in  circular  orbits  about  their  common  center  of  mass,  that  their  angular 
frequency  (0)  is  given  by  ft2  **  G(M^  +  M2) /a3,  and  that  the  P  particle 
remains  in  the  plane  in  which  and  M2  revolve.  In  the  rotating  refer¬ 
ence  frame  whose  origin  is  the  center  of  mass  and  whose  angular  frequency 
is  ft,  the  two  components  of  the  equation  of  motion  of  a  g  particle  with 
coordinates  (X,Y)  are 
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ABSTRACT 

The  classical  restricted  three-body  problem  is  generalized  to  include  the  force 
of  radiation  pressure  and  the  Poynting-Robertson  effect.  The  positions  of  the 
Lagrangian  points  and  are  found  as  functions  of  &,  the  ratio  of  radiational  to  ' 
gravitational  forces.  The  Poynting-Robertson  effect  renders  the  and  L,.  points 
unstable  on  a  time  scale  (T)  long  compared  to  the  period  of  rotation  of  the  two 
massive  bodies.  For  the  solar  system,  T  is  given  by  T  «=  [ (l-pJ^/B)  544  a2  yr  where 
a  is  the  separation  between  the  Sun  and  the  planet  in  AU.  Implications  for  space 
colonization  and  a  mechanism  for  producing  azimuthal  asymmetries  in  the  interplanetary 
dust  complex  are  discussed;  the  latter  can  be  observationally  tested  by  the  Zodiacal 
Light  Experiment  aboard  the  International  Solar  Polar  Mission  spacecraft  to  be  launched 
in  1983. 


